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Abstract 

Using the Jordan- Wigner transformation and the continued-fraction method we calculate exactly the 
density of states and thermodynamic quantities of the periodic nonuniform spin-i isotropic XY chain in a 
transverse field. We discuss in detail the changes in the behaviour of the thermodynamic quantities caused 
by regular nonuniformity. The exact consideration of thermodynamics is extended including a random 
Lorentzian transverse field. The presented results are used to study the Peierls instability in a quantum spin 
chain. In particular, we examine the influence of a non-random/random field on the spin-Peierls instability 
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1 Introduction 



The study of regularly nonuniform spin models is an attracting problem of statistical mechanics. Besides of 
its general academic importance the development of magnetic materials in recent years makes the study of 
nonuniform spin models particularly interesting for experimental application. In order to achieve a progress in 
understanding the generic features generated by periodic nonuniformity it is desirable to examine simple models 
that may be investigated without making any approximation. Among possible candidates of such systems one 
can mention the spin-i XY model in one dimension. The uniform one-dimensional spin-i XY model in a 
transverse field was introduced by Lieb, Schultz and Mattis. 1 These authors noticed that a lot of statistical 
mechanics calculations for such a spin model can be performed exactly since it can be rewritten as a model 
of noninteracting spinless fermions by means of the Jordan- Wigner transformation. The nonuniform version of 
the transverse spin-i XY chain also can be mapped onto a chain of free spinless fermions, however, with an 
on-site energy and hopping integrals that vary from site to site. Especially attractive is the case of isotropic 
spin coupling with regularly alternating exchange integrals and transverse fields since after fermionization one 
faces a model for which a lot of work has been done. One should mention here the results for the tight-binding 
Hamiltonian of periodically modulated chains 2 ' 3 and the spinless Falicov-Kimball chain. 4,5 

One of the goals of the present paper is to give a magnetic interpretation of those results derived for the 
one-dimensional tight-binding spinless fermions. Exploiting the continued-fraction approach developed in the 
mentioned papers 2-5 we shall be able to calculate exactly the one-fermion Green functions and therefore to 
obtain the thermodynamic quantities for the periodic nonuniform spin-i isotropic XY chain in a transverse 
field. We shall treat few examples of the periodic nonuniform spin-i isotropic XY chain in a transverse field 
in order to reveal the changes in the thermodynamic properties induced by periodic nonuniformity. The model 
of the considered regularly nonuniform magnetic chain allows even a natural extension to include additional 
disorder remaining the model exactly solvable. Namely, one can assume the transverse fields to be random 
independent Lorentzian variables with regularly alternating mean values and widths of distribution. To derive 
exactly the random-averaged density of states for such a model one should at first average a set of equations 
for the Green functions using contour integrals. 6-11 As a result one comes to a set of equations similar to that 
for the periodic nonuniform non-random case. 

It should be noted that the periodic nonuniform spin-i isotropic XY chain was considered in several 
papers 12-25 dealing mainly with the adiabatic treatment of the spin-Peierls instability. However, those papers 
were focussed mostly on the influence of the structural degrees of freedom upon the magnetic ones, rather than 
on the exhaustive analysis of the properties of a magnetic chain with regularly alternating exchange couplings. 
Another closely related study concerns the spin-i isotropic XY model on a one-dimensional superlattice 26 . The 
treatment reported in Ref. 26, however, was restricted to the magnon spectrum. A related study of the spin-| 
isotropic XY chain in a transverse field with two kinds of coupling constant aimed on examining the condition 
for appearance of an energy gap was reported in Ref. 27. Finally, let us note that the periodic nonuniform chain 
can be viewed as the uniform chain with a crystalline unit cell containing several sites of the initial lattice (as a 
matter of fact such a point of view was adopted, for example, in Refs. 12, 13) and thus the standard methods 
elaborated for such complex crystals may be exploited. However, we prefer to treat periodically nonuniform 
chains since from such a viewpoint an elegant continued-fraction approach immediately arises that seems to be 
a natural and convenient language for describing such compounds. 

The present paper is a more extensive version of the results briefly reported in Refs. 28, 29 containing more 
details on the calculation and more applications. We show that the presented below method based on continued- 
fraction representation for the one-fermion diagonal Green functions immediately reproduces the results for the 
spin-i transverse isotropic XY chains having periods 2 and 3 and in contrast to other approaches easily yields 
the results for larger periods (e.g., 4 and 12). It should be stressed that the elaborated approach is a systematic 
method that permits to consider in the same fashion the regularly nonuniform chains with randomness that 
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cannot be done within the frames of the approaches exploited in Refs. 12-27. We present a comprehensive 
study of the thermodynamic properties (density of states, gap in the energy spectrum, entropy, specific heat, 
magnetization, susceptibility) of regularly alternating chains having periods 2, 3, 4, and 12 discussing in detail 
the dependences of the energy gap and the low-temperature transverse magnetization on the transverse field 
and comparing the latter dependence with the corresponding one for the classical spin chain. We underline a 
possibility of a nonzero transverse magnetization at the zero average transverse field in the periodic nonuniform 
chain owing to a regular nonuniformity. This material constitutes Section 2. In Section 3 we demonstrate the 
influence of a diagonal disorder on the effects caused by periodic nonuniformity assuming the transverse fields 
to be independent random Lorentzian variables. For simplicity we restrict ourselves by the case of a chain 
having period 2. The results obtained for the spin-i transverse isotropic XY chains with period 2 are applied 
for an analysis of the spin-Peierls instability in the adiabatic limit with respect to dimerization in the presence 
of a non-random/random (Lorentzian) transverse field (Section 4). We find how the (random) transverse field 
influences the dependence of the (averaged) total energy on the dimerization parameter tracing a suppression 
of dimerization by the non-random (random) field. 



2 Periodic nonuniform spin-^ isotropic XY chain in a transverse 
field 

Let us consider a cyclic nonuniform isotropic XY chain of N (eventually N — > oo) spins s = | in a transverse 
field. The Hamiltonian of the system reads 

N N 

h = J2 a,** + 2 ]>> + a yy n+1 ) 

n—l n—1 

N / i\ N 

= fin ( S n S n ~ ^ ) + 51 ^" ( S n S n+l + S n S n+l) > S n+N = S n~ (1) 

71=1 ^ ' 71=1 

Here tt n is the transverse field at site n and 2I n is the exchange coupling between the sites n and n + 1. Let 
us note that s z = J2n=i s n commutes with the Hamiltonian H (1) and hence the eigenfunctions of H can be 
classified according to eigenvalues of s z . Moreover, at 0„ = the ground state of H corresponds to s z = 0. 
After making use of the Jordan- Wigner transformation one comes to a cyclic chain of spinless fermions governed 
by the Hamiltonian 

N ( 1\ N 

H = ^2n n ( c+c„ - 2 ) + H 7 " ( c « c «+i _ c « c rt+i) • ( 2 ) 

n=l ^ ' n=\ 

The so-called boundary term is not essential for calculation of thermodynamic functions 30 and has been omitted. 
We shall discuss the most general case, i.e. assuming that both transverse fields and exchange couplings vary 
from site to site. Note, that in the particular case when the transverse field is uniform one recognizes in Eq. 
(2) the Hamiltonian of the system considered in Ref. 2. In addition, in another limiting case after substitution 
Cl n — > Uw n , I n — > —t Eq. (2) transforms into the Hamiltonian of a one-dimensional spinless Falicov-Kimball 
model in the notations used in Refs. 4, 5. 

Let us introduce the temperature double-time Green functions 

Gnm(t) = Ti0(±t)({c n (t),c+(O)}), Gl m {t) = (l/27r)/_ 00 oo dccxp(-i^)GT m ( w ±k), e - +0, where the 
angular brackets denote the thermodynamic average. Consider further the set of equations of motion for 
G7un = GJ m (w±ie) 

(u±ie- fl n ) G^ m - I n - 1 G n z _ lm - I n G^ +lrn = 5 nm . (3) 
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Our task is to evaluate the diagonal Green functions GJ„, the imaginary part of which gives the density of 
states p(u)), 



1 N 



(4) 



that on its part, yields the thermodynamic properties of the spin model (1). 

It is a simple matter to obtain from Eq. (3) the following representation for GJ n 



nn 



l 



lo ± ie - fi„ - A„ - A£ ' 

J n-1 



w ± ie — 



aj±le — f2 n _2- 



A+ 



uj ± ie — O n+ i 



aj±le — On+2 



(5) 



Equations (4), (5) are extremely useful for examining thermodynamic properties of the periodic nonuniform 
spin-^ isotropic XY chain in a transverse field, since the evaluation of periodic continued fractions 31 emerging 
in (5) is quite simple and reduces to solving quadratic equations. 

It should be noted here that the continued-fraction representation of the one-particle Green functions has 
been widely used for tight-binding electrons over the last two decades. As an example let us refer here to 
the papers of Haydock, Heine and Kelly 32,33 and the review articles. 34,35 However, those studies were aimed 
mainly on getting the electronic band structure of non-translationally invariant systems (alternatively to the 
band theory) starting from the local environment of atom and in practice were connected with an appropriate 
approximative termination of continued fractions. In what follows we shall use the exact values of continued 
fractions (since they are periodic) to reveal the effects of regular nommiformity on the magnon band structure. 

Consider at first a uniform chain QqISIqI .... In this case one comes to a periodic continued fraction having 
a period 1 



A+ 



lu ± ie — Qq 



i 2 



1 ± ie - Q - A ' 



(6) 



The quadratic equation for A (6) can be solved with 
A = 



± ie - n + \]{lo±i6- n f - AP 
and therefore p(u) according to (4), (5) becomes 

p(uj) -- 



0, 



■ ± ie - fio - \/(w±ie-fto) 2 -4/ 2 j 



(7) 



if 4/ 2 - (w - n ) 2 > 0, 

otherwise. 



(8) 



The self-consistent equation for the continued fraction (6) introduces a spurious root. However, the false solution 
is eliminated requiring p(w) to be not negative. Let us emphasize the attractive features of the continued-fraction 
approach reminding how p(u>) (8) can be obtained within the frames of the standard technique. Usually one 
substitutes into Eq. (3) GJ m = (1/N) ^ K exp[i(n — to)k]GJ to obtain GJ = 1/(uj ± ie — Oo — 27 cos k) and 
then evaluates the integral GJ„ = (l/2n) da / (u ± ie — fio — 2/cosk) using, for example, contour integrals 

to get GJ n = l/J(u) ± ie — £! ) 2 — 4/ 2 and therefore the density of states (8). 
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The advantages of the continued-fraction approach become clear while treating the periodic nonuniform 
chains. We shall demonstrate this in some detail for regularly modulated chains with periods of modulation of 
2, 3 and 4. 

(i) Consider a regular alternating chain Clili^h^ih^h • • • • In this case periodic continued fractions 
having a period 2 emerge. Solving similar quadratic equations as (6) for A~, A+, A~ +1 , A+ +1 one obtains as 
a result the Green functions GJ„, G^ +1 n+1 and therefore the density of states p(u) 

„„,_< A^^. if %)>»■ 

0, otherwise; 

1 2 



B{u) = - [(w - fii) (w - o 2 ) - A 2 - /|] ' 

= -(u- h) (lu - b 2 ) (w - 63) (w - 64) . (9) 



Here 61 > b 2 > b 3 > 64 denote the four roots of the equation B(to) = 0, namely 



\- 2 <-l. ■ il 2 : I b,. i (Ox + fi 2 ) ± b 2 l r (10) 



withbi = ^ (Qi - n 2 ) 2 + 4 (\h\ + \I 2 \)\ b 2 = £y(fii -fi 2 ) + 4 (|/i| -|/ 2 |) . Solving the inequality > 
one can write the density of states (9) in the explicit form 

if LU < 64, 63 < lu < b 2l bi < U!, 



27r v 8 M 

The result for the uniform chain (8) is contained in the density of states (11), (10), (9) as a partial case when 
fii = 2 = too, h = h = I- 

(ii) Next we consider the regularly modulated chain 01/102/2^3-^3^1^1^2-^2^3-^3 ■ • ■ • In this case one 
generates the corresponding periodic continued fractions of period 3. Going along the lines as described above 
one gets 

0, otherwise; 
C(uj) = 4i?if if - [if (lu - 3 ) + II (w - fii) + / 2 (w - fi 2 ) - (w - Oi) (w - 2 ) (w - 3 )] 2 

6 

= n ( W Cj ) , (12) 

where Cj are the six roots of the equation C(uu) = 0. To find them one must solve two qubic equations that 
follow from Eq. (12). 

(iii) Finally, let us consider the regularly modulated chain O1/1O2/2O3/3O4/4O1/1O2/2O3/3O4/4 .... In this 
case one gets periodic continued fractions with period 4. The density of states for such a chain is given by 

r ^^m, if v(u)>o, 

p(cu)=\ 47r V^) 

I 0, otherwise; 

W(w) = l\ (2lo - 3 - H 4 ) + J 2 2 ( 2w - fii - H 4 ) + / 3 2 (2w - Oi - 2 ) + if (2w - 2 - 3 ) 

— (w — Qi)(w — 2 )(w — 3 ) — (lu — Oi)(w — 2 )(w — O4) 

— (lu — Oi)(w — 3 )(w — 4 ) — (w — Q 2 )(w — 3 )(w — O4), 

V(lu) = 4/i 2 / 2 2 / 3 2 /l - [(w - Oi) (w - 2 ) (w - 3 ) (w - 4 ) 

-/1 2 (lu - 3 ) (w - 4 ) - / 2 2 (w - Oi) (lu - 4 ) 



■5 



-jf (w - Oi) ( w - n 2 ) - / 4 2 (w - n 2 ) (w - n 3 ) 
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+ ^7| + /|/|] 2 = -J]( w - dj .), 



(13) 



where dj are the eight roots of the equation T>(ui) = 0. To find them one must solve two equations of 4th 
order that follow from Eq. (13). Let us note that all dj (as well as all Cj) are real since they can be viewed as 
eigenvalues of symmetric matrices. 2 

There are no principal difficulties in proceeding the analytic calculations of p(tu) for larger periods, except 
the fact that they become more cumbersome. All Green functions required for getting the density of states p(u>) 
(4) are calculated by solving quadratic equations, however, further analysis of the band structure is becoming 
more complicated. This analysis, however, can be easily implemented on a computer and the results for the 
chains having period 12 presented below were obtained in such a manner. 

Let us discuss the results for the density of states for the considered periodic nonuniform chains. The main 
consequence of introducing the nonuniformity is a splitting of the initial magnon band into several subbands 
(compare (8) and (9) - (13)). The edges of the subbands are determined by the roots of equations B(w) = 0, 
C(ui) = 0, T>(u>) = 0, etc.. p(uj) is positive inside the subbands, tends to infinity inversely proportionally to the 
square root of u) — ui e when u approaches the subbands edges u e , and is equal to zero outside the subbands. The 
number of subbands does not exceed the period of the chain. At special (symmetric) values of the Hamiltonian 
parameters the roots of the equation that determines the subband edges may become multiple and the zeros in 
the denominator and the numerator in the expression for p(u>) may cancel each other. As a result due to an 
increase of symmetry one may observe a smaller number of subbands. This 'mechanism' is easily traced, for 
example, in formulas (9) - (11) if putting fix — Q 2 , \Ii\ — \I 2 \- The described magnon band structure can be 
seen in Figs. 1 and 2a, 3 and 4a, and 5 and 6a where we show p(u>) for a few particular periodic nonuniform 
chains having periods 2, 3, and 12, respectively. The splitting caused by periodic nonuniformity in fact is not 
surprising. The periodic nonuniform chain is simply another viewpoint on the uniform chain with a crystalline 
unit cell containing several sites. On the other hand, it is generally known that one may expect several subbands 
for a crystal having several atoms per unit cell. 36 

Further one can easily calculate the widths of energy gaps in the magnon spectrum that appear due to 
nonuniformity. For example, for a chain having a period 2 one finds 



This quantity is connected with a gap between the ground state energy and the first excited state energy of the 
spin chain. As an example consider the chain Oo/i51o^2^o-^i^o^2 • • • • The edges of the upper magnon subband 
are given by flo + |7i| + \I 2 \ and flo + \ \Ii \ — |7 2 ||, whereas the edges of the lower magnon subband are given 
by flo — I \h I — I-Z2 1 1, and flo — \h\ — l-^l- At fl — the ground state (for which s z = 0) corresponds to the 
filled lower subband and the empty upper subband and the energy spectrum exhibits a gap A(0) = ||7i| — I/2II 
(this is the energy required to create a hole in the lower subband — the first excited state of the spin chain 
(with s z 7^ 0)). With increasing of flo the gap A decreases as A(f2o) = ||/i| — IJ2II — ^0 and becomes zero at 
flo = I \Ii I — l-^H- With further increasing of flo the gap remains equal to zero up to the value of the transverse 
field flo = \h\ + \h\ after which the gap opens and increases as A(f2 ) = flo — \h\ — \h\ (the ground state 
for the transverse field larger than |ii| + \I 2 \ corresponds to the empty subbands and the written A(f2 ) is the 
energy required to create a particle in the vicinity of the lower edge of the lower subband). For chains with 
larger periods one finds more complicated behaviour of the energy gap A with varying of the field flo (see Fig. 
7a and Figs. 7b - 7d). 

The splitting of the magnon band into subbands caused by nonuniformity has interesting consequences for 
thermodynamic properties. The entropy, specific heat, transverse magnetization and static transverse linear 



b 2 -b 3 = v /(fi 1 -0 2 ) 2 +4(|7 1 |-|7 2 |) 2 . 



(14) 
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susceptibility are determined through the density of states according to the following formulas 



dEp(E) 



, E \ E E 
In 2 cosh — — — tanh ■ 



2kT J 2kT 2kT 



E 
2kT 



(15) 



m 



Xz 



— OO N / 

OO 

> = d£p(£)tanh^, (17) 

— OO 
OO 

-iff / d ^>I^i7X- < 18 > 



Apparently the most spectacular changes caused by regular nonuniformity are observed in the dependence 
of transverse magnetization (17) on transverse field at low temperatures (Figs. 2b, 4b, 6b). Since for T —* 0, 
tanh 2^ tends either to — 1 if E < 0, or to 1 if E > 0, one immediately finds due to the splitting of the 
magnon band into subbands that the low-temperature dependence of m z versus Oo must be composed of 
sharply increasing parts (they appear when E = moves with increasing of O from the bottom to the top 
of each subband) separated by horizontal parts (they appear when E = moves with increasing of O inside 
the gaps). The number of plateaus is determined by the number of subbands. It should be emphasized here 
that a study of magnetization plateaus for quantum spin chains is a hot topic at the present time. 37 However 
in such studies usually more general spin chains are attacked which cannot be treated within the frames of the 
described approach. For example, the spin-i XXX chain can be mapped onto the chain of interacting spinless 
fermions with the intersite interaction of the same order as the hopping integral and hence the results derived 
rigorously for noninteracting fermions cannot be immediately extended for this more complicated spin chain. 

It is interesting to note that the appearance of plateaus in the dependence of transverse magnetization on 
transverse field at T = for the regularly nonuniform isotropic XY chains essentially differs in the quantum 
and classical cases. The Hamiltonian of the classical nonuniform isotropic XY chain in a transverse field reads 

N N 

cos9 n + 2 Ins 2 cos(0„ — <pn+i) sin sin#„ + i (19) 

n—l n—1 

that immediately yields the ansatz for the ground state energy in the uniform case 

E = Nn a scos6- 2N\I\s 2 sm 2 6 = N(l Q m z + 2N\I\ (m 2 z - s 2 ) (20) 

where the ground state transverse magnetization m z = s cos 9 has been introduced. Minimizing E with respect 
to cos# one finds that for s = | the quantity — m z increases as while f2 increases from to 2\I\ and 

— m z = i with further increase of fig. Using numerical calculations for finite chains (the number of spins TV is a 
multiple of 12) with periods 2, 3, 12 we found that the detailed profiles for the quantum and classical chains are 
different, although the values of the transverse field at which a saturation of the transverse magnetization occurs 
are the same. Though one could argue that the magnetization plateaus are connected with the quantum nature 
of the spins we found for special parameter sets even in the classical chain plateaus in the dependence —m z 
versus Oo (compare dashed curves in Figs. 8a, 8b and in Figs. 2b, 4b). For instance, the well pronounced plateau 
shown by the dashed line in Fig. 8b occurs at the same height as in the quantum case. The corresponding 
classical state is a state ITTITT4TTITT-- - where the arrows symbolize classical spins pointing either in —z- or 
+z-dircction. An evident difference between the quantum and classical case is connected with the slope of the 
m z (O ) curve at T — 0. The slope remains finite in the classical case but becomes infinite approaching the 
plateaus in the quantum case. The infinite slope in the quantum case is clearly a consequence of the singularities 
in the density of states. 
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One of the interesting magnetic properties of the periodic nonuniform spin-^ isotropic XY chain is the possi- 
bility of the existence of a non-zero transverse magnetization m z at zero average transverse field (X)^=i ^« = 0)- 
For illustration we consider as an example a chain having the period 4 and the parameters Oi = Q3 = 0, 
0,2 = — £1 4 < 0, |ii| = |/2| > 0, I/3I = I/4 1 = 0. At site n + 1 we have the transverse field U 2 < surrounded on 
the left and right side by the strong couplings |7i| = \I2\- At site n + 3 we have the transverse field — f2 2 > 
surrounded by the weak couplings I/3I = I/4I = 0. One may expect that the local transverse magnetization 
at site n + 1 has a smaller value and opposite direction with respect to that quantity at site n + 3 and there- 
fore a non-zero total transverse magnetization at zero average transverse field may be expected. Consider the 
described chain in more detail. From Eq. (13) for the above set of parameters it follows that 

p(u) = \i5\u! v — + X 2 S (u>) + X36 lu> z— + A 4 d (u + 2 ) (21) 



and the coefficients Xj may be found comparing (21) and (13) in the vicinity of ^zVgi±g£i , , [Wgs+gfi and 
—0,2- As a result one gets Xj = j (see Fig. 9a). Now transverse magnetization (17) at T = is m z = — | 7^ 
although — (solid curves in Figs. 9b, 9c; in the latter figure the solid curve especially clearly shows 

that — m z = g). If |/ 3 | = \I 4 \ ^ the magnon subbands look as in Fig. 9a and at T = one has m z — (Figs. 
9b, 9c). However, such a position of the subbands provides an interesting temperature dependence of m z at 
J2n=i — (dashed and dotted curves in Fig. 9c) reminding the 'order from disorder' phenomenon, 38-40 i.e. 
increasing of order with increasing temperature. 

Let us turn to other thermodynamic quantities. Every infinite slope in the dependence m z versus f2o at 
T = induces a singularity in the dependence Xzz versus fto at T = 0. However, there is no need to plot this 
dependence. Since l/4fcTcosh 2 ^ tends to 5(E) as T -> one gets from (18) that at T = -\zz = p(0). 
The latter dependence as a matter of fact can be seen in Figs. 2a, 4a, 6a. The changes in the temperature 
dependences of entropy and specific heat due to nonuniformity which are displayed in Figs. 2c, 4c, 6c and 2d, 
4d, 6d can be understood while bearing in mind the behaviour of integrands in (15), (16) that are products of 
the functions with evident dependences on the temperature and the density of states. Note that as a result of 
the magnon band splitting the temperature dependence of the specific heat may exhibit a two-peak structure 
(Fig. 2d) or even a more complicated behaviour (solid curve in Fig. 4d). Finally we look at Xzz- As mentioned 
above at T = we have — Xzz — p(0). Analysing the density of states depicted in Figs. 2a, 4a, 6a one 
finds that nonuniformity may either suppress or enhance the initial (that is at Qq = 0) static transverse linear 
susceptibility — \zz at T — shown in Figs. 2f, 4f, 6f. 



3 Periodic nonuniform spin-^ isotropic XY chain in a random Lo- 
rentzian transverse field 

In this Section we consider a generalization of model (1) including additional randomness in the transverse 
fields. We assume the transverse fields to be independent random variables each with a Lorentzian probability 
distribution 

P( Q n) = l- £V-=i- ' (22) 

7T(Q „-Q n ) +T 2 n 

Here f2o« is the mean value of the transverse field at site n and T n is the width of its distribution. We are 



interested in the random-averaged density of states p(u>) that follows from the random-averaged diagonal Green 
functions G% n according to Eq. (4). Repeating the arguments presented in Refs. 6-11 one gets the following 
set of equations for the random-averaged Green functions 



(w ± ir„ - fi 0n ) Gnm - I n -lGn-l,m ~ ^n^n+l,m — ^nm (23) 
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that immediately yields 



1 n-l 



W i ir„_i — Qo.n-l 



w±ir„_ 2 -o 



2 — ^0,71-2" 



A+ = ^ js • (24) 

w ± iT n+ i - n , n+ i - ^r^T^— 

In case fion, r n , J„ vary regularly from site to site one again comes to the periodic continued fractions. They 
can be calculated as solutions of the corresponding quadratic equations. Thus one gets rigorously the random- 
averaged Green functions and therefore the random-averaged density of states. For example, for a regular 
random chain Q iTiIiil 2^2l2^oi^iIi^02^2l2 . . . one finds 

yn = (r, + r 2 )^M+ZH _ sgnB>)(2w _ n 01 - n oa) ^- g/ ^ , 

= vW)) 2 + (S"(c)) 2 , 

B'(u>) = [(w - o i)(w - r> 02 ) - rxr 2 - / 2 - /|] 2 - [(w - Q i)r 2 + (w - f> 02 )ri] 2 - 4/ 2 / 2 , 

= 2 [(w - noi)(w - ft 02 ) - - / 2 - 1$] [(to - n 01 )r 2 + (w - o 02 )ri] . (25) 

The random-averaged density of states (25) transforms into (9) if ri = T 2 = 0, and into the result reported in 
Ref. 9, rfu) = T(l/7r)Iml/v/(o) ± ir - fi ) 2 - 4/ 2 , if Q i = «o 2 = fi , I\ = T 2 = T, h = h = I. 

Let us discuss the effects of the considered diagonal Lorentzian disorder. The main effect of the randomness 
is smearing out the band structure. However, one can see a difference in smoothed magnon subbands for the 
uniform disorder (when Ti = r 2 ) (see Fig. 10a) and the nonuniform disorder (when T± =/= r 2 ) (see Fig. 11a). 
Namely, in the former case both subbands are smeared out in the same way, whereas in the latter case, the 
subbands are smeared out differently and, at least for small strengths of disorder, in one subband the peaks at 
the band edges persist. This circumstance in the latter case induces an interesting step- like behaviour of the 
low-temperature transverse magnetization as a function of transverse field. Namely, as can be seen in Fig. lib 
the disorder smooths only one step in contrast to Fig. 10b in which both steps are smeared out. The difference 
in the influence of the uniform and nonuniform disorders on other thermodynamic quantities can be seen in 
Figs. 10c - lOf and 11c - llf. 



4 Periodic nonuniform spin-^ isotropic XY chains and spin-Peierls 
instability 

In this Section we want to demonstrate that the results for the density of states of the periodic nonuniform 
spin-i isotropic XY chains obtained within the continued-fraction approach may be of use for the study of the 
spin-Pcicrls instability in these chains in adiabatic limit. The discovery of existence of the spin-Peierls transition 
in the inorganic compound CuGeOs ' has stimulated much research work in this field. In particular, the 
influence of an external field or randomness attracts much interest both from experimental and theoretical 
viewpoints (see e.g. Refs. 42-49). 

Let us start from the non-random case. In order to examine the instability of the spin chain with re- 
spect to dimerization one must calculate the ground state energy per site of the regularly alternating chain 
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^1/1^2/2^1/1^2/2 ■ • ■ (sec Eqs. (9) - (11)) 

e = --j dE P (E)\E\ = -- j dE-= 



-b 2 



(\E- 



fti + \E + n\ 



V- (£ 2 -b 2 )(£ 2 -b 2 ) 

where ft = (fti + ft2)/2. Depending on the value of ft formula (26) can be rewritten as follows 

-b 2 



(26) 



eo 



dE 



\n\\E\ 



yj- (E 2 -ti{)(E 2 -b 2 2 ) 



(27) 



if bi < |ft|, 



eo 



-|fi| 

"i / « 



£ 2 



-b 2 



V-(/; 2 -b 2 )(^-b|) 



d/; 



Iftll^l 



-|n| 



V- (£ 2 - b?) (£ 2 - b 2 ) 



(28) 



if b 2 < |ft| < bi, and 



eo 



4/ 



£ 2 



V- (£ 2 - b 2 ) (£ 2 - b 2 ) 



(29) 



if |ft| < b 2 . Introducing a new variable ip by the relation E 
final expression for the ground state energy 

1 

71 



-\Jb\ — (b\ — b 2 ) sin 2 tp one gets the following 



eo 



b 2 



b? 



Ift 



(I 



biE ^, 

where E(?/>,a 2 ) = dipy/l — a 2 sin 2 93 is the elliptic integral of the second kind 50 and 

0, 



(30) 



arcsm 



b 2 -n 2 



if bi < |ft|, 
if b 2 < |ft| < bi, 
if |ft| < b 2 . 



(31) 



The result obtained by Pincus 13 follows from (30), (31) if fti = ft2 = 0. However, the described approach 
permits to get the ground state energy (or the Helmholtz free energy) for more complicated regular nonuni- 
formities (e.g., for chains with regularly alternating non-random or random (Lorcntzian) transverse fields). 
To demonstrate this let us consider at first the spin-Peierls instability with respect to dimerization in the 
presence of a non-random transverse field. We introduce dimerization parameter S and assume in (30), (31) 
|/i| = |/|(1 + S), IZ2I = |/|(1 — 5), < S < 1. Taking into account that the elastic energy per site is aS 2 one 
must seek the minimum of the total energy £(S) = e Q (S) + aS 2 as a function of 5. For £ (6) we find 



£{5) = 



V(r!i - fi 2 ) 2 + 16/ 2 



2tt 



EU, T 



4/ 2 (l - S 2 



(fti -ft 2 ) 2 + 4/ 2 



-|fti+ft 2 



2ir 



aS 2 



(32) 



with 



0, 



arCSm Y 4/ 2 (l-<5 2 ) i 



2 ' 



if v/(fti - ft 2 ) 2 + 16/ 2 < |fti + ft 2 |, 

if y/tfh. - ft 2 ) 2 + 16/ 2 <5 2 < |fti + ft 2 | < v/(fti - ft 2 ) 2 + 16/ 2 , 
if |fti + ft 2 | < y/(Qi - ft 2 ) 2 + 16P5 2 . 



(33) 



Eqs. (32), (33) in the limit of uniform field fti = ft 2 coincide with the result reported in Ref. 16. For strong 
fields |fti + ft 2 | > y/Wi - ft 2 ) 2 + 16/ 2 one finds that £{$) = -±|fti + ft 2 | + aS 2 and the equation = 
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has only the zero solution 5* = (no dimerization in strong enough fields), whereas for weaker fields besides 
the zero solution there may be a non-zero one 5* ^ coming from the equation 



= y/jsi! - n 2 y + 16/ 2 

a ~ 4tt(1 - S 2 ) 



4/ 2 (l-<5 2 ) \ ( AI 2 {\-5 



F ^'i(fl 1 - n 2 ) 2 + 4/2 J E - fi 2 ) 2 + 4/2 



(34) 



where F(-0, a 2 ) = dip/ \f\ — a 2 sin 2 ip is the elliptic integral of the first kind. 50 

In the following discussion of results we choose a uniform transverse field Q\ = Q 2 = f2o, |^o| < 2|/|. To 
give a guide for further reading this paragraph we summarize the main results valid for sufficiently hard lattices 
(having a > ^-). (i) For zero field we have a minimum of the total energy £(S) at a nonzero value of the 
dimerization parameter S* ^ 0. (ii) For finite but small fields £(S) still exhibits one minimum at 5* ^ the 
position of which remains unchanged, (iii) When the field achieves a certain characteristic value flo a a second 
local minimum appears at 6* = 0. The two minima at 5* = and S* are separated by a maximum, (iv) At 
a second characteristic field fiofc both minima at S* = and S* ^ have the same depth, (v) Further increasing 
fio the minimum at 5* = becomes the global one and at a certain characteristic field £Iq c the minimum at 
S* 7^ abruptly disappears. The scenario described in (i) - (v) is typical for a first order transition characterized 
by the order parameter 5* and driven by the transverse field Oo- Now we illustrate it in a more detail. 

In Fig. 12 we show for different values of a how the dependence of £ (5) — £(0) on the dimerization parameter 
varies with the strength of the field CIq. As it follows from Eqs. (32), (33) (and can be also seen in Fig. 12 
where, however, the difference £{8) — £(0) is depicted) the total energy £{8) at sufficiently large values of S 
(S > ^jj, £(S) — £(0) at the value is denoted by dark circles in Fig. 12) becomes independent of the field. 
In Fig. 13 we plot the solution of Eqs. (34), (33) for different lattices (i.e. different values of a) in the presence 
of the field. As a matter of fact we calculated r.h.s. of Eq. (34) varying 8 from to 1 and finding in such a way 
for what a this value of 6* realizes. Note that solutions of Eqs. (34), (33) S* which are smaller than realize 
a maximum of the total energy, whereas solutions 5* which are larger than realize a minimum. This can 
be seen, for example, for a lattice with a = 0.4 in Figs. 12c and 13b, 13c: at f^o = 0-1 the total energy £(§) 
exhibits two minima at 5* = and 5* =/= separated by a maximum at intermediate value of 5*; at fig = 0.2 
the total energy £(S) exhibits only a minimum at S* = 0. From Figs. 12, 13 and Eqs. (34), (33) one concludes 
that for soft lattices having a < ^ there is no solution of Eqs. (34), (33) fulfilling the presupposition S* < 1. 
Such lattices are excluded from further consideration. For other lattices the solution of Eqs. (34), (33) 5* ^ 
existing for zero transverse field does not feel the presence of a small field, however, abruptly vanishes at a 
certain value of the transverse field. Moreover, for soft lattices one needs larger fields than for hard lattices for 
a disappearance of the solution of Eqs. (34), (33) (compare Figs. 13b - 13f with Fig. 13a). Thus, in the case 
of hard lattices even small transverse fields may destroy the dimerization. As it is seen e.g. for a lattice with 
a = 0.2 (Figs. 12, 13) above a certain characteristic value of the transverse field f2oa (for which Eqs. (34), 
(33) has the solution 8* = 0) (Ooo ~ 0.2) £{S) starts to exhibit in addition to the global minimum at S* ^ 0, 
a local one at S* = 0, two minima are separated by a maximum at the intermediate value of the dimerization 
parameter. With increasing of Qo the depths of the minima at first become equal (when Qo na s a characteristic 

value flob) and then the minima at 5* — becomes a global one. The latter minima remains the only one at 

inol 



Slo having a characteristic value flo c (for which Eqs. (34), (33) has the solution 5* = i^jj) (Qq c w 0.5) that 
manifests a complete suppression of the dimerization by the field. In Fig. 14 we show different regions in the 
plane transverse field Oo - lattice parameter a in which £ (5), < 6 < 1 exhibits one minimum at 5* ^ 
(region A), two minima at 6* = and 5* =/= separated by a maximum (regions Bi and B2; in the region Bi 
the minimum at S* =/= is deeper, whereas in the region B2 the minimum at S* = is deeper), one minimum at 
5* = (region C). To find the line that separates Bi and B 2 one must find for a given O such a 6* at which 
£{8) — £(0) (32), (33) with a given by the r.h.s. of Eq. (34), (33) equals to zero, and then to evaluate the r.h.s. 
of Eq. (34) at the sought S* . Crossing the phase diagram by a vertical line corresponding to a certain lattice 
(e.g. with a = 0.2 in Fig. 14) one obtains the field at which the first order transition between the dimerized 



fl 



and uniform phases occur (Slob in Fig. 14) and the width of hysteresis (determined by S7oa and Oo c in Fig. 14). 

Next we consider the influence of a random Lorentzian transverse field on the spin-Peierls instability with 
respect to dimerization. For that we calculate the difference in random-averaged total energy (to avoid non- 
physical infinities due to the Lorentzian probability distribution) 

oo 

W)-W) = -l J dE(^(E)-ME))\E\+a8 2 (35) 



with ps(E) given by Eq. (26) where |ii| = |/|(1 + 5), I/2I = — &)■ Let us start from the case Sloi = ^02 = 0, 
Ti = T2 = r generalizing in such a way the consideration for the zero transverse field by assuming the latter 
to be random (Lorentzian) with the zero mean value. As can be seen in Fig. 15 the randomness leads to a 
continuous decrease of the non-zero value of dimerization parameter at which the random-averaged total energy 
exhibits minimum. At sufficiently large strengths of disorder T the minimum of the random-averaged total 
energy occurs already at the zero dimerization parameter, i.e. randomness acts against dimerization and may 
suppress it completely for sufficiently large strength of disorder. Considering the equation 



9£(S) _ 1 f ^dps(E) 
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\ J ™ d -^\E\ + 2aS = (36) 



one can find its solution 8* for different T (see Fig. 16). From Fig. 16 one sees that in the case of hard lattices 
even small disorder may destroy the dimerization. In Fig. 17 we depicted different regions in the plane strength 
of disorder T - lattice parameter a in which £(8) — £(0), < 8 < 1 exhibits one minimum at 8* 7^ (region 
A) or one minimum at 8* = (region C). The boundary curve between the regions C and A is obtained by 
calculating a from (36) with varying T for fixed 6 = 0. Thus, the random field with zero mean value suppresses 
dimerization with increasing the strength of disorder, however the dimerization parameter 8* vanishes according 
to a second order phase transition scenario in contrast to the previous case. 

Finally we consider the case of random field with non-zero average value, i.e., Sloi = ^02 = !!o / 0, 
Ti = T 2 = T. For small strengths of randomness T the above discussed scenario of one or two minimum in 
£(S) in dependence of the value of the field remains valid. A switching on randomness for a system being in the 
region A at T = (Fig. 14) leads to continuous decreasing of 8* ^ to zero. For a system being in the regions 
Bi or B2 an increasing of randomness usually leads at first to a continuous decrease of 8* 7^ with a decrease 
of the depth of that minimum and then to an abrupt disappearance of 8* 7^ above a certain strength of 
disorder. We also observed another influence of small randomness for a system being in the region Bi, namely, 
an increasing of randomness leads at first to a disappearance of the minimum at 8* = that appears again for 
larger strength of disorder. The details can be traced in Fig. 18 where we plotted the dependence £(8) — £ (0) vs 
8 for different T considering two mean values of the random transverse field Slo = 0-1 and ^0 = 0.3 and in Fig. 
19 where wc illustrated the vanishing and appearance of the minimum at 8* = with increase of randomness. 
Both the one minimum profile (solid curve in Fig. 18b) and the two minima profile (solid curves in Figs. 18c, 
18e) of that dependence existing in the non-random case T = are finally destroyed by increasing disorder. The 
phase diagrams in the T - a plane for the two mentioned values of Slo are shown in Fig. 20. 

Closing this Section, we want to make some comments concerning the conclusions on spin-Peierls instability 
that can be drawn using exact results for thermodynamic quantities of regularly nonuniform spin-i isotropic 
XY chain in a transverse field. Although the described basic picture of a first order phase transition in a uniform 
field seems to be qualitatively correct we should keep in mind that an increasing of field at low temperature leads 
to a transition from dimerized to incommensurate phase. This fact was observed experimentally and analysed 
theoretically mainly for the models of CuGeOa in a number of papers. 42 ' 51-54 Clearly, the simple ansatz for 
the lattice distortion S1S2S1S2 ■ ■ ■ , Si + 82 = permitted us to compare the ground state energies only for 
dimerized and uniform phases. To detect a transition from the dimerized to the incommensurate phase with 
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increasing of field one may analyse the ground state energy of a chain having larger period, say 12. The presence 
of randomness requires even more complicated lattice distortions to be examined and the continued-fraction 
approach for rigorous study of thermodynamics of the regularly alternating spin-i isotropic XY chain in a 
transverse field provides some possibilities to perform such an analysis. We must also keep in mind that the 
known spin-Peierls compounds are described by the spin-i isotropic Heisenberg chain rather than XY chain, 
however, one may expect that the basic features of the studied phenomenon should be similar for both quantum 
spin models. 

5 Summary 

To summarize, we have studied rigorously the magnon density of states and the thermodynamics of the periodic 
nonuniform spin-^ isotropic XY chain in non-random/random (Lorentzian) transverse field. We have exploited 
the Jordan- Wigner transformation, the temperature double-time Green functions and the continued fractions. 
The Green functions approach seems to be the most convenient tool for a study of thermodynamics of the 
considered spin chains since it permits to examine such models with regular nonuniformity or some type of 
randomness or both. Regular nonuniformity leads to a splitting of the magnon band into subbands that 
in its turn leads to some spectacular changes in the behaviour of the gap in the energy spectrum and the 
thermodynamic quantities. In particular, the low-temperature dependence of the transverse magnetization on 
the transverse field is composed of sharply increasing parts separated by plateaus, the temperature dependence 
of specific heat may exhibit a well pronounced two-peak structure, the temperature dependence of the initial 
transverse linear susceptibility may be enhanced or suppressed. Regularly nonuniform spin-^ isotropic XY 
chain may exhibit a non-zero transverse magnetization at the zero average transverse field. The regularly 
alternating Lorentzian disorder in the transverse field may in specific manner influence the thermodynamic 
quantities leading, for instance, to a smearing out of only one 'step' in the step-like dependence of the transverse 
magnetization versus the transverse field at T = 0. The derived results for the (random-averaged) ground state 
energy permit to analyse the effects of external non-random/random field on the spin-Peierls instability. Both, 
magnetic field as well as randomness may destroy the dimcrization as the analysis of the (random-averaged) 
total energy manifests. 

The presented treatment of the regularly periodic spin-i isotropic XY chains is restricted to the density of 
states and therefore only to thermodynamics. It will be interesting to study the effects of periodic nonuniformity 
on spin correlations and their dynamics especially for a model of spin-Peierls instability. Some work for the 
dynamic zz spin correlations for such models has been done in Ref. 16. Another interesting problem concerns 
the treatment of the periodic nonuniform spin-i transverse XY chains with an anisotropic exchange coupling 
(and in particular the extremely anisotropic case, i.e. the spin-i transverse Ising chain). Some results for 
thermodynamics of such regularly nonuniform chains having period 2 were obtained in Refs. 14, 17, 23. Their 
relation to the spin-Peierls instability seems to be an intriguing issue. 
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List of figure captions 



FIG. 1. Magnon band structure for periodic chains f2i/if2 2 / 2 f2i/if2 2 /2 • • ■ , = Slo + ty; the shadowed 
areas correspond to the allowed magnon energies, a) + Vt' 2 = 2, \Ii\ = \I 2 \ = 0.5; b) tt[ + Sl 2 = 2, \Ii \ = 0.75, 
\I 2 \ = 0.25; c) n[ = Q 2 = 1, \h\ + \I 2 \ = 1; d) Q[ = 1.5, Q 2 = 0.5, |/i| + \I 2 \ = 1. The horizontal lines single out 
the following particular chains: il[ = fi 2 = 1, |ii| = \I 2 \ = 0.5 (dotted curves), f2'i = 2, fl 2 = 0, |/i| = \I 2 \ = 0.5 
(dashed curve), fl[ = Ct' 2 = 1, |/i| = 0.75, \I 2 \ = 0.25 (dashed-dotted curves), Q[ = 1.5, n' 2 = 0.5, |ii| = 0.75, 
\I 2 \ = 0.25 (solid curves). 



FIG. 2. The density of states (a), the dependence of the transverse magnetization on transverse field at 
T = (b), the temperature dependence of the entropy (c), specific heat (d), transverse magnetization (e), and 
static linear transverse susceptibility (f) at £1 = for periodic chains f2i/ifi 2 / 2 f2i/iSl 2 / 2 ■ • ■ 7 % = ^0 + Q'j- 
The dotted curves correspond to the uniform case Cl[ = Cl 2 = 1, |ii| = I/2I = 0.5, the dashed curves correspond 
to the case Cl[ = 2, Q! 2 = 0, |7i| = I/2I = 0.5, the dashed-dotted curves correspond to the case fi'i = Q! 2 = 1, 
|/i| = 0.75, \I 2 \ = 0.25, and the solid curves correspond to the case £l[ = 1.5, Q,' 2 = 0.5, |ii| = 0.75, \I 2 \ = 0.25. 



FIG. 3. Magnon band structure for periodic chains f2i/i£l 2 / 2 f2 3 / 3 £li/iil 2 / 2 £l 3 / 3 • • ■ , = ^0 + f^-; the 
shadowed areas correspond to the allowed magnon energies, a) Oj + fl 2 + Q 3 = 3, Q[ — 2fl 2 + Q 3 = 0, 
l^i I = l-k I = \h\ = 0.5; b) fl^+Qi. + n's = 3, Q[-2n' 2 + Q 3 = 1.5, |/i| = |/ 2 | = |/ 3 | = 0.5; c) Qi+ft 2 + Q 3 = 3, 
Q[ - 2Q' 2 + Q' 3 = 0, |/i| = 0.75, \I 2 \ = 0.5, |/ 3 | = 0.25; d) Q[ + W 2 + Q' 3 = 3, Q[ - 2ft 2 + Q' 3 = 1.5, |/i| = 0.75, 
\I 2 \ - 0.5, |/ 3 | = 0.25; e) fl[ = Q 2 = Q 3 = 1, \h\ + \I 2 \ + \h\ = 1.5, \h\ - 2\I 2 \ + 1/ 3 | = 0; f) O'j = 2 = Q 3 = 1, 
\h\ + \h\ + \h\ = 1.5, |/i| - 2|/ 2 | + |/ 3 | = 0.75; g) n[ = 1.5, 2 = 1, 3 = 0.5, \h\ + \I 2 \ + \I 3 \ = 1.5, 
|/i| - 2|/ 2 | + |/ 3 | = 0; h) fl[ = 1.5, n> 2 = 1, 3 = 0.5, |/i| + |/ 2 | + |/ 3 | = 1.5, |/i| - 2|/ 2 | + |/ 3 | = 0.75. The 
horizontal lines single out the following particular chains: Cl^ = Cl 2 = Q 3 = 1, |/i| = |/ 2 | = |/ 3 | = 0.5 (dotted 
curves), n[ = 2.5, fi 2 = 0.5, fi 3 = 0, |/i| = |/ 2 | = |/ 3 | = 0.5 (dashed curve), fl[ = 1, 2 = 0.5, fl' 3 = 1.5, 
|/i| = 0.75, \I 2 \ = 0.5, \I 2 \ = 0.25 (dashed-dotted curve), Q[ = 1.5, 2 = 1, 3 = 0.5, |/i| = 1, \I 2 \ = |/ 3 | = 0.25 
(solid curve). 



FIG. 4. The same as in Fig. 2 for periodic chains f2i/i>T2 2 / 2 f2 3 / 3 £li/if2 2 / 2 £l 3 / 3 • • • , = ^0 + ty- The 
dotted, dashed, dashed-dotted, and solid curves correspond to the cases pointed out in the capture to Fig. 3. 



FIG. 5. The same as in Figs. 1, 3 for periodic chains having a period 12, flih . . . ^12/12^1/1 • • • ^12*12 ■ • ■ , 
fli = fl 2 = . . . = Qq, f2 7 = fig = . . . = fii2, I\ = I 2 = . . . = Ie, I7 = *8 = • • • = Ii2, = ^0 + fij-- a) 
n[ + Q' 7 = 2, |/i| = |/ 7 | = 0.5; b) fl[ + Q' 7 = 2, |/i| = 0.75, |/ 7 | = 0.25; c) Q[ = Q' r = 1, |/i| + |/ 7 | = 1; d) 
f2'i = 1.5, S1' 7 = 0.5, |/i| + I/7I = 1. The horizontal lines single out the following particular chains: fl[ = Q! 7 = 1, 
|/i| = |/ 7 | = 0.5 (dotted curves), fl[ = 2, fl' 7 = 0, |/i| = |/ 7 | = 0.5 (dashed curve), Q[ = Q' 7 = 1, |/i| = 0.75, 
|/ 7 | = 0.25 (dashed-dotted curves), n[ = 1.5, Q' 7 = 0.5, |/i| = 0.75, \I 7 \ = 0.25 (solid curves). 
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FIG. 6. The same as in Figs. 2, 4 for the chains singled out in Fig. 5. 



FIG. 7. The dependence of the energy gap A between the ground state and the first excited state on 
transverse field fio for certain regularly nonuniform chains, a) The chain fio-Zifio/ 2 fio7ifio/ 2 • • • > \Ii\ = 0.75, 
|J 2 | = 0.25; b) - d) the chains having periods 2, 3, and 12, respectively, with the notations as in Figs. 2, 4, 6. 



FIG. 8. The dependence of the transverse magnetization on the transverse field fi at T — for classical 
periodic nonuniform isotropic XY chains in a transverse field, a) Chains having a period 2 (fii = 2, Vl' 2 = 0, 
l^i I = \h\ = 0.5 (dashed curve), Q[ = Q' 2 = 1, \h\ = 0.75, \I 2 \ = 0.25 (dashed-dotted curve), Q[ = 1.5, 
fi 2 = 0.5, |7i| = 0.75, \I 2 \ = 0.25 (solid curve)); b) chains having a period 3 (fii = 2.5, fi 2 = 0.5, fi!, = 0, 
\h\ = \h\ = I/3 1 = 0.5 (dashed curve), Q[ = 1, fi 2 = 0.5, fi 3 = 1.5, |Ji| = 0.75, \I 2 \ = 0.5, |J 3 | = 0.25 (dashed- 
dotted curve), fii = 1.5, fi 2 = 1, fi 3 = 0.5, |7i| = 1, \I 2 \ = \h\= 0.25 (solid curve)); c) chains having a period 
12 (fii = 2, fi' 7 = 0, |ii| = \h\= 0.5 (dashed curve), fi' x = fi' 7 = 1, = 0.75, \h\ = 0.25 (dashed-dotted 
curve), fii = 1.5, fi 7 = 0.5, = 0.75, \h\ = 0.25 (solid curve)). 



FIG. 9. Illustration of the existence of a non-zero transverse magnetization at the zero average transverse 
field in a chain having period 4. fii = fi^ = 0, fi^ = — fi^ = — 1, |ii| = II2I = 0.5, I/3I = I/4I = (solid curves), 
|/ 3 | = \h\ = 0.05 (dashed curves), |7 3 | = \h\ = 0.25 (dotted curves). 



FIG. 10. The random-averaged density of states (a), the dependence of the transverse magnetization on 
transverse field at T = (b), the temperature dependence of the entropy (c), specific heat (d), transverse magne- 
tization (e), and static linear transverse susceptibility (f) at fio = for periodic chains 

fioiri/ifio2r 2 / 2 fioir 1 / 1 fio 2 r 2 / 2 ..., fi 0j = fi + fii = 1.5, fi 2 = 0.5, = 0.75, |/ 2 | - 0.25 for the 

case of uniform disorder ri = T 2 = T. The solid curves correspond to the non-random case r = 0; the 
long-dashed curves correspond to T = 0.1; the short-dashed curves correspond to T = 0.25; the dotted curves 
correspond to T — 0.5. 



FIG. 11. The same as in Fig. 10 for nonuniform disorder Ti ^ 0, T 2 = 0. The solid curves correspond to the 
non-random case Ti = 0; the long-dashed curves correspond to Ti = 0.1; the short-dashed curves correspond 
to Ti = 0.25; the dotted curves correspond to T\ = 0.5. 



FIG. 12. Change of the total energy £(S) — £(0) as a function of the dimerization parameter S in the presence 
of the uniform transverse field; |7| = 0.5; a) a = 0, b) a = 0.2, c) a = 0.4; fi = (solid curves), fio = 0.1 
(dashed-dotted-dotted curves), fio = 0.2 (dashed-dotted curves), fio = 0.3 (dashed curves), fi = 0.4 (dotted 
curves) . 
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FIG. 13. Dimerization parameter 5* as a function of a in the presence of a uniform transverse field Qq; 
\I\ = 0.5; Q a = (a), fl = 0.1 (b), Q a = 0.2 (c), fl = 0.3 (d), fl = 0.4 (e), fl = 0.5 (f). The solid curves 
show the solution of Eqs. (34), (33) corresponding to a minimum of the total energy; the dashed curve in (a) 
corresponds to the dependence S* versus a valid for hard lattices that was obtained in Ref. 13; the dashed 
curves in (b) - (f) show the solution of Eqs. (34), (33) corresponding to a maximum of the total energy. 



FIG. 14. Different types of solution for the dimerization parameter 5* (0 < 5* < 1) in the plane f2 - 
a; \I\ = 0.5. Region A: £(5) has one minimum at S* ^ 0, regions Bi, B2: £{$) has two minima at 5* = 
(favourable in B2) and 5* 7^ (favourable in Bi) separated by a maximum, moreover, the depths of the minima 
at the line that separates Bi and B2 are the same; region C: £(5) has one minimum at 5* = 0. 



FIG. 15. Change of the random-averaged total energy as a function of the dimerization parameter in the 
presence of a uniform random Lorentzian transverse field with zero mean value; |/| = 0.5, Ti = T2 = T = 
(solid curves), T = 0.02 (dashed-dotted curves), T = 0.1 (dashed curves), T — 0.5 (dotted curves); a) a = 0, b) 
a = 0.2, c) a = 0.4. 



FIG. 16. The solution of Eq. (36) as a function of a in the presence of disorder; |/| = 0.5, O01 = ^02 = 0, 
ri = T 2 = T = (solid curves), T = 0.02 (dashed-dotted curves), T = 0.1 (dashed curves), T = 0.5 (dotted 
curves). 



FIG. 17. Different types of solution for the dimerization parameter S* in the plane T - a; |/| = 0.5, Oo = 0. 
Region A: £(S) — £(0) has one minimum at S* ^ 0, region C: £{5) — £(0) has one minimum at S* = 0. 



FIG. 18. Change of the random-averaged total energy as a function of the dimerization parameter in the 
presence of the uniform random Lorentzian transverse field with a non-zero mean value Oo = 0.1 (a, b, c) and 
ft = 0.3 (d, e, f); |7| = 0.5, T 1 = T 2 = T = (solid curves), T = 0.02 (dashed-dotted curves), T = 0.1 (dashed 
curves), T — 0.5 (dotted curves); a = (a, d), a = 0.2 (b, e), a = 0.4 (c, f). 



FIG. 19. Change of £(S) — £(0) as a function of 5 in the presence of the uniform random Lorentzian transverse 
field with flo = 0.3, T = 0.01 (solid curves), T = 0.1 (dashed-dotted curves), T = 0.2 (dashed curves), T = 0.3 
(dotted curves); |/| = 0.5, a = 0.15. 



FIG. 20. Different types of solution for the dimerization parameter 5* in the plane in the plane T - a; 
|/| = 0.5, n = 0.1 (a), n = 0.3 (b). Region A: £jS) - £j0) has one minimum at 5* ^ 0, region Bi: £jS) -£{0) 
has two minima at 5 ^ and 5 = and the first one is favourable, region B 2 : £ (S) — £(0) has two minima at 
5^0 and 5 = and the second one is favourable, region C: £(5) — £(0) has one minimum at 5* = 0. 
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Figure 1: FIGURE 1. 
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Figure 2: FIGURE 2. 
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Figure 4: FIGURE 4. 
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Figure 5: FIGURE 5. 
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Figure 6: FIGURE 6. 
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Figure 7: FIGURE 7. 
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Figure 8: FIGURE 8. 
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Figure 9: FIGURE 9. 
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Figure 10: FIGURE 10. 



28 




Figure 11: FIGURE 11. 
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Figure 12: FIGURE 12. 
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Figure 13: FIGURE 13. 
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Figure 14: FIGURE 14. 
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Figure 15: FIGURE 15. 
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Figure 16: FIGURE 16. 
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Figure 17: FIGURE 17. 
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Figure 18: FIGURE 18. 
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Figure 19: FIGURE 19. 
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Figure 20: FIGURE 20. 
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